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We show that an extension of the standard BCS Hamiltonian leads to an infinite number of 
condensates with different energy gaps and self-similar properties, described by a cyclic RG flow of 
the BCS coupling constant which returns to its original value after a finite RG time. 

PACS numbers: ll.10.Hi, 74.20.Fg, 71. 10. Li 



The Renormalization Group (RG) continues to be one 
of the most important tools for studying the qualitative 
and quantitative properties of quantum field theories and 
many-body problems in Condensed Matter physics. The 
emphasis so far has been mainly on flows toward fixed 
points in the UV or IR. Recently, an entirely novel kind 
of RG flow has been discovered in a number of systems 
wherein the RG exhibits a cyclic behavior: after a finite 
RG transformation the couplings return to their original 
values and the cycle repeats itself. Thus if one decreases 
the size of the system by a specific factor that depends 
on the coupling constants, one recovers the initial sys- 
tem, much like a Russian doll, or quantum version of the 
Mandelbrot set. Bedaque, Hammer and Van Kolck ob- 
served this behavior in a 3-body hamiltonian of interest 
in nuclear physics . This motivated Glazek and Wilson 
to define a very simple quantum-mechanical hamiltonian 
with similar properties 0. In the meantime such be- 
havior was proposed for a certain regime of anisotropic 
current-current interactions in 2 dimensional quantum 
field theory @. 

The models in jl], |2| are problems in zero-dimensional 
quantum mechanics, and are thus considerably simpler 
than the quantum field theory in || . In the latter, stan- 
dard quantum field theory methods of the renormaliza- 
tion group were used, however knowledge of the beta 
function to all orders was necessary to observe the cyclic 
flow. What is somewhat surprising is that the model 
considered in H is not very exotic, and is in fact a well- 
known theory that arises in many physical problems: at 
one-loop it is nothing more than the famous Kosterlitz- 
Thouless RG flow, where the cyclic regime corresponds 
to \g±\ > \g»\. This motivated us to find a simpler 
many-body problem that captures the essential features 
of the cyclic RG behavior. We found that a simple exten- 
sion of the BCS hamiltonian has the desirable properties. 
Namely, our model is based on the BCS hamiltonian with 
scattering potential Vjj/ equal to g + i9 for Ej > £y and 
g — i6 for Sj < sy in units of the energy spacing 8. 

The main features of the spectrum are the following. 
For large system size, there are an infinite number of 
BCS condensates, each characterized by an energy gap 
A n which depends on g, 9. The role of these many con- 
densates becomes clearer when we investigate the RG 
properties. As in the models considered in pi 0, ffl, the 



RG flow possesses jumps from g = +oo to g = — oo and 
a new cycle begins. Let L = e~ s Lo denote the RG scale, 
which in our problem corresponds to N the number of 
unperturbed energy levels, and A the period of an RG 
cycle: g(e~ x L) = g(L). We show that A = n/9. 

The model we shall consider is an extension of the 
reduced BCS model used to describe ultrasmall super- 
conducting grains although our results are valid for 
more general cases. Let c t ± (cj.±) denote creation- 
annihilation operators for electrons in time reversal states 
|±). The index j = 1,..,N refers to N equally spaced 
energy levels Sj with — ui < ej < u. The energy sj rep- 
resents the energy of a pair of electrons in a given level. 
The level spacing will be denoted 25, i.e. £j+i — £j = 25, 
so that uj — N5 is twice the Debye energy. The ener- 
gies in this paper are twice their standard values pi . Let 
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I) j = Cj t -Cj t +, M = c'j + c'j _ denote the usual Cooper- 
pair operators. The Hilbert space TLn is spanned by the 
combination of empty and occupied states. At half-filling 
the dimension of the Hilbert space is the combinatorial 
number C^ 2 . 

Our model is defined by the reduced BCS hamiltonian 

JV N 

where Vjy is the scattering potential. In the usual BCS 
model, Vjj' is taken to be a constant. Here we add an 
imaginary part which breaks time reversal, 



Vjj> = < G 



G + i@ if Sj > £f 

if £j = Ey 



G-iQ 
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G = gS 
Q = 9S 
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This hamiltonian is hermitian since V*y = Vyj. We con- 
sider the positive dimensionless couplings g and 9. 
The BCS variational ansatz for this model is 
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The mean-field treatment yields the well-known equa- 
tions 



Ei = 



= I e 2^ (i-k. 

Cj' = £j — j-! — Vjj , 



(4) 
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where Aj and <fij satisfy the gap equation: 



A, 
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The chemical potential equation in our case is satisfied 
with fj, = 0. In the thermodynamic limit, N — > oo and 
5^0, with fixed o> = 7V<5, the sums over £j become the 
integrals de/28. The gap equation turns into 



r de' A(e') 
A(e) = o / 7-4 +t6 
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where A(e) = A(e)e^ £ ). 

Differentiating (||) with respect to e yields 

d<j) e 



de E(e) ' 
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and the condition that A(e) = A is independent of e. 
The solution to eq. ([?]) can be taken to be 



0(e) = flsinh -1 — . 
Using eq. (§h in the gap equation (pi) gives 



(8) 



f >(u ' d(b 9 
1= — (gcoscj) + Osincj)) => tan0(w) = -. (9) 
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Solving eq. (0) for the gap yields an infinite number of 
solutions A„. They can be parameterized as follows: 

u) nix 
A « = -T7-' t n = t + —, n = 0,1,2,..., (10) 
sinn t n 9 

where to is the principal solution to the equation 



Next we derive RG equations for our model. Let c/jv, 
On denote the couplings for the hamiltonian Hn with 
N energy levels. The idea behind the RG method is 
to derive an effective hamiltonian if/v-i depending on 
renormalized couplings <?jy— i, Qn-i by integrating out 
the highest energy levels En or e\. This can be accom- 
plished by a canonical transformation, which is formally 
analogous to the one used to derive the t — J model from 
the Hubbard model at strong coupling 0. 

We perform the calculation for general V. The inte- 
gration of the level e n yields 



V n> ~ V n> + 2 V V Nj , t- , 

(13) 

where £j = Sj — Vjj. Integration of the level e\ gives the 
same eq. ( fL3[ ) with the replacement £jv — £,j —> — £i + £j- 
Specializing to the potential eq. (|^) and approximating 
£N~£j or — Si+Sj by u — NS, the above equation implies 



AN) 



9N-1 = fjN + Jf(a% + °n)> 



7JV-1 — 0]y. 



(14) 



Thus 9 is unrenormalized. 

In the large N limit one can define a variable s = 
log No /N, where Nq is the initial size of the system. Then 
the beta function reads 

The solution to the above equation is 

ff ( S ) = 0tan[0 S + tan- 1 (^)] , g = g(N ). (16) 
The main features of this RG flow are the cyclicity 



tan(^o) 



< to < 



29' 



(11) 



The gaps satisfy A > Ai > • • •. Each gap A„ repre- 
sents a different BCS eigenstate |v4cs)- O ne can show 

that |(4osl4cs)l < ex P [~ N ( A n - A n ,) 2 /8u 2 ] , in the 
limit where A„ <C u>. Thus in the large N limit, these 
eigenstates are orthogonal and should all appear in the 
spectrum, together with the usual quasi-particle excita- 
tions above them. In the limit 9^0 the gaps A„>o — > 0, 
and since to — l/g, Ao ~ 2we -1 / 9 , in the weak coupling 
regime, recovering the standard BCS result. 

For weak coupling models A„ <C w, all the gaps are 
related by a scale transformation A„ ~ 2N5e~ to ~ n7T / e . 
Therefore, defining the condensation energy of the n-th 
BCS eigenstate as E^ = (V'bcsI^IV'bcs) — ^fs we get 



E. 



(n) 
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Thus the spectrum of condensation energies reflects the 
scaling behavior of the gaps. 



g(s + X) =g(s) 



g(e- x N)=g(N), A 



(17) 



and the jumps from +oo to — oo, when reducing the size. 

The cyclicity of the RG has some important implica- 
tions for the spectrum. Let {E(g,9, N)} denote the en- 
ergy spectrum of the hamiltonian Hn- The RG analysis 
implies we can compute this spectrum using the hamil- 
tonian H]y/(g(N')) if g(N') is related to g(N) according 
the RG equation (|1|). Moreover, if N' and N are re- 
lated by one RG cycle, N' = e~ x N, then g(N') = g(N). 
Thus a plot of the spectrum {E(g,9, N)} as a function 
of N but at fixed g, 9 is expected to reveal the cyclicity 
{E(g, 9, e- x N)} = {E(g, 9, AT)}. Since our RG procedure 
is not exact, we expect to observe this signature within 
the range of our approximations, i.e. for \E\ <C u). In- 
deed, this agrees with the result shown in eq. (fT2|). This 
can also be observed in fig. 1 for the one Cooper pair 
case, with the cyclicity given by Ai = 2A (see below). 

Eliminating go in eq. (|l^) in terms of the mean-field 
solution, eqs. ( |Toj . [TT| ) , we observe that the jumps in g(s) 
from +00 to -co occur at scales s = t n . As N decreases, 
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g increases steadily to +00 and then jumps to — 00. At 
g = +00, to — 0, t\ = Tt/0,..., whereas for g = —00, 

t Q = tt/9, ti = 2n/9, Plugging this into eq. (|l0|), one 

readily sees that 



A (.9 



-00) 
-00) 



00 

A n ( 5 



-00), 



(18) 



which indicates that at every jump the lowest condensate 
disappears from the spectrum, since Eq ' (g = +00) = 
—00. Eq. (fL8j) implies, for the remaining condensates, 



E { £ +X \g = +00) = E^{g = -00). This result is 
in agreement with eq. (|l2|). Therefore, the condensate 

I^BCS ) °f one cycle plays the same role as V'bcs) 
of the next cycle. 

The blow up of Ao and Eq at g = +00 is an artifact 
of the RG scheme used here, since we can only trust 
the RG for energies below the cutoff u>. However the 
disappearance of bound states is correctly described by 
this RG (see the one Cooper pair problem for a more 
detailed discussion). 

When N = 00 the infinite number of condensates are 
all expected to appear in the spectrum. However at finite 
TV this is not possible since the Hilbert space is finite 
dimensional. One can use the RG to estimate the number 
of condensates nc in the spectrum as a function of N. 
From the discussion above a condensate disappears from 
the spectrum for each RG cycle. Thus nc should simply 
correspond to the number of cycles in log N: 



nc 



log AT. 



(19) 



So far we have found a close relationship between the 
spectrum of our extended BCS hamiltonian in the mean- 
field approximation and the RG flow of the coupling con- 
stants. In order to get a further confirmation of our re- 
sults we should compute the spectrum for a finite size 
system. However, it is very difficult to reach intermedi- 
ate sizes for this model numerically, since the dimension 
of the Hilbert space grows as 2 N /N 1 / 2 . Fortunately, to 
this end, the similarities between the many-body case 
and the case of one Copper pair, in the presence of the 
Fermi sea, are widely known. 

For one Cooper pair in the presence of the Fermi sea, 
consider an eigenstate of the form \ip) — J2j V'j&jlO)- The 
Schrodinger equation reads 

fa - E)^ = GiPj + {G + iQ) ^ ^ + (G - i6) 

Kj l>j 

(20) 

with £j £ (0,w), i.e. the Fermi sea is not accessible for 
the pair. Then, in the large TV limit, the sums ^ - are 
replaced by integrals L de/25, leading to 



dp' 

(s E)^e) =g I — ^e') + 1 



(21) 



Differentiating the above with respect to e and inte- 
grating, one obtains 



1 



Jdlog(e-E) 



e-E 



(22) 



This wave function does not have cuts in two cases, i) 
E < and ii) E > 10, Case i) corresponds to the usual 
Copper pair problem where one is looking for bound state 
solutions (these are the solutions we claim to have a sim- 
ilar behavior to the many-body case). Plugging eq. ( p^ ) 
back into eq. (j2l|) one finds 



9 + it 
9 -it 



= 1- 



E 



(23) 



This equation has an infinite number of solutions given 
by 



E n — — - 



tn — to 



2nn 



n e 2, 



e*» - 1 : 

where to is the principal solution to the equation 
'1 



tan I -#<oJ = -, 



< t Q < 



(24) 



(25) 



The n > solutions correspond to E n < 0, while those 
with n < yield E n > lo. 

As for the many-body case the spectrum has a scaling 
behavior for weak coupling systems, namely 



E n 



-N5 e - t o- 2m / , 



(26) 



An RG analysis similar to the one in the many-body 
problem leads to the equation 



9n-2 — 9n + 



9 2 N + e\ 



N 



>N-l 



(27) 



N-l-gtr' 
which in the large N limit becomes 

g = i^ + 0. (28) 

The solution to this equation is given by eq. ( |l6| ) just by 
replacing 9s — > Os/2. This implies that the period of the 
cyclicity in s — XogNo/N is Ai = 2tt/9. 

The factor 1/2 in the above formulas, as compared to 
the many-body problem, comes from the non accessibility 
of the Cooper pair to the states below the Fermi level. 

The discussion leading to eq. (|l^) can be repeated for 
the one-Cooper pair problem where the role of the n-th 
condensate is played by the n-th bound state with energy 
E n given in eq. (24), obtaining 



E (9 

E n +i(g 



-00) = —00 

-00) = E n (g = -00). 



(29) 



Thus we expect that in each RG cycle a bound state 
will disappear. The analogue of equation ( |l9|) is 



n B 



9 , N 

2^ l0g Y' 



(30) 



4 



where ns is the number of bound states in the spectrum. 

This again shows the agreement between the mean- 
field and the RG results. We confirm bellow this picture 
with numerical calculations. 




FIG. 1: Exact eigenstates of one-Cooper pair Hamiltonian 
for N levels, from No = 500 down to 50. We depict only 
the states nearest to zero. The vertical lines are at the values 
N n = e~ nAl No- The dotted horizontal lines show the cyclicity 
of the spectrum. 



Fig. 1 shows the numerical solution of eq. (|2C|) for g = 
1. 9 = 12 and N ranging from 500 down to 50. For each 
N there are ns{N) bound states E n < 0, where Ub(N) 
is in good agreement with eq. (p0|). 

The spectrum shows the self-similarity found in the 
approaches above: scaling the system by a factor e~ Al , 
with Ai = 2it/9, one recovers the same spectrum for suf- 
ficiently small energies, i.e. 



E n+1 (N,g,9) =E n (e- x ^N,g,9) 



(31) 



Fig. 1 also shows the existence of critical values N cn , 
in the intervals (e~ nAl N, e~(" +1 ) Al N), where the bound 
state closest to the Fermi level disappears into the "con- 
tinuum" . This effect leads to the reshuffling of bound 
states, n + 1 — > n, observed in eq. (|3l|). The critical 
sizes are also related by scaling, i.e. N c _ n /N c _ n+ i = e~ Xl . 
All these phenomena are in good agreement with the RG 
interpretation we proposed, where the condensates dis- 
appear at scales where g(s = t n ) = Too. All these points 
are related by the scaling factor e~ Xl . 

The RG behavior is presented in fig. 2, which shows the 
eigenvalues E n (N) of the one-Cooper pair Hamiltonian, 
with <7tv running under eq. (p7|). The spectrum remains 
unchanged for E n (N) <C NS, as shown in fig. 2b. In 
fig. 2a one observes that for the energies E n (N) > NS 
the result of the RG is not reliable. Nevertheless the 
RG flow describes qualitatively the disappearance of the 
lowest bound state and the reshuffling of energy levels 
after a cycle, and furthermore at the predicted scales. 



In summary, we have shown that adding to the stan- 
dard BCS Hamiltonian a time reversal breaking term, 
parametrized by a coupling constant 9, generates an infi- 
nite number of condensates with energy gaps A„ related, 
for weak BCS couplings g, by a scale factor e~ A with 
A = tt/9. This unusual spectrum is explained by the 
cyclic behavior of the RG flow of g, which reproduces 
itself after a finite RG time s equal to A. We have also 
solved the finite temperature BCS gap equation, obtain- 
ing a critical temperature T Cj „ for the n-th condensate 
which is related to the zero temperature gap A n (0) ex- 
actly as in the BCS theory, i.e. A n (0)/T Ci „ = 3.52 for 
weak couplings. 

The simplicity of the model proposed in this letter sug- 
gests that Russian doll superconductors could perhaps be 
realized experimentally. Finally, we point out that the 
critical temperature can be raised by varying 9. 



to 
g 
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FIG. 2: Eigenstates of one-Cooper pair Hamiltonian with <?jv 
given by eq.(^) with go = 1 and 9 = 12. The vertical lines 
denote the positions at which g jumps from +oo to — oo. 
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